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Abstract 

We prove the existence of non-decaying real solutions of the Johnson equation, van- 
ishing as x — * +oo. We obtain asymptotic formulas as t — > oo for the solutions in the 
form of an infinite series of asymptotic solitons with curved lines of constant phase 
and varying amplitude and width. 

1 Introduction 

The Johnson equation (JE) 

1 3 v \ Via 2 

V t + ~V XXX + ~VV X + = ^-Uyy (1.1) 

(a 2 = ±1) or the cylindrical Kadomtsev-Petviashvili equation, is the analogue of the well- 
known cylindrical Korteweg-de Vries equation (a = 0) in two spatial dimensions (2D). 
The JE was obtained firstly in M under the description of the surface waves on a shallow 
incompressible liquid. Later it was shown that it describes the propagation of waves in 



the stratified media pi. It follows from the derivation of (1.1) that the correct statement 
of the Cauchy problem is possible only as t = to > 0. 

In H and Q the equivalence of the Kadomtsev-Petviashvili equation (KP) and the JE 
was established. Let u(£, r\, r) be an arbitrary solution of the KP 

1 3 \ 3a 2 , . 

Ur + + -jUU/: j = ^-u m . (1.2) 

Then the function 

v(x,y,t) =u(x- j,t\ (1.3) 
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satisfies the JE. This mapping 77, r) — ► v(x,y,t) is invertible. Each solution v(x,y,t) 
of the JE generates a solution of the KP by the formula 

»«, ,,r)=«(«+ 3^, fA (1.4) 



It was shown in y] that mappings (|l.3| ) and (|l.4j ) preserve the class of functions rapidly 
decaying at infinity (as (x 2 + y 2 )~ 1 — ► 0), and all the results obtained in the theory of the 
KP in the corresponding class of solutions can be directly applied to solve the Johnson 
equation. The situation is similar for the solutions involving the Airy functions (see || and 
||). Obviously this is not the fact for periodic initial data, which are not invariant with 
respect to this transformation, and investigation of the JE is an independent interesting 



problem in this case (for the KP see the corresponding theory, for example, in [35|-p 



We are interested in the construction of a class of JE-I (a = i in (1.1)) non-decaying 
solutions, which are bounded for all (x, y, t) and vanish as x — > +oo for all fixed y and 
t. Such a kind of KP solutions was constructed and investigated firstly for KP-II in M, 



HI, and then for KP-I in |J-[13|. It turns out that all basic stages of the construction of 
the solutions of the KP and the JE, and the study of their asymptotic behaviour admit 
mutual recounting using the described mapping. 
We apply the change of variables 

y 2 t yt 



48a 2 ' V 4 ' 



to the scheme of the V.E.Zakharov and A.B.Shabat "dressing method" [14] of integration 
of the KP, and obtain analogous formulas for the JE. Using them we prove the existence 
of a class of JE-I non-decaying solutions with the prescribed properties. The simplest one 
is the one-soliton solution 

t 2q2 
v(x,y,t) 



cosh 2 



(p S R, q G R ) which corresponds to the KP plane-soliton by virtue of (|1.3|). 

Then we study the asymptotic behaviour of the constructed solution as t — > oo. The 
investigations of long-time asymptotic behaviour of non-decaying solutions of 2D non- 
linear evolution equations is closely connected with the same investigations in one spatial 
dimension. A.V.Gurevich and L.P.Pitaevsky studied in 1973 a non-decaying solution of 
the Korteweg-de Vries (KdV) equation, which describes the evolution of an initial step- 
function ( (lf| , [l(| ) . They applied Whitham method to construct an approximation of this 
solution by a knoidal wave with slowly varying parameters and detected the appearance 
of many strong oscillations like solitons on the front of the solution for a large time. This 
approximate solution satisfies the KdV-equation with error vanishing as t — ► +oo. The 



mathematical ground of this phenomenon was done in 1975 by E.Ya.Khruslov in [17] and 



[18], where the nature of these solitons was explained. Subsequently these solitons were 
called asymptotic solitons. An analogous phenomenon of splitting of non-decaying initial 
data into infinite series of solitons was proved later for other KdV-like equations (nonlinear 
Schrodinger equation, sine-Gordon equation, modified KdV and the Toda lattice as a 
discrete analogue of the KdV) 
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In @-fLf] the method proposed by E.Ya.Khruslov (@,|8|) was extended to the in- 



vestigation of the asymptotic behaviour of non-decaying solutions of KP-type equations 
(KP, modified KP-I and 2D-Gardner equation) as t — > oo. It was proved that they are 
represented as infinite series of solitons with curved lines of constant phase in the neigh- 
bourhood of the front as t — > oo. These asymptotic solitons were called curved asymptotic 
solitons. Note that recently (f30|) V.E.Zakharov also considered a curved soliton of the 
KP-II equation, but in another space-time domain. 

Our principal goal is to prove the phenomenon of splitting of non-decaying solutions of 
the JE-I into infinite series of solitons as t — > oo. 

We prove that there exist non-decaying real solutions of the JE-I, which split in the 
neighbourhood of the front into a series of solitons of the form 

v n (x 1 y,t) = — T ? ^ ■ (1.5) 



cosh 



q (y) (x - C(y)t + ^ (lni»+V2 _ l ng ( y ) _ n ( y ))) 



which depend on two parameters C(y) and g(y). The functions po(y), qo(y) and 4> n (y) are 

n + l/2 

completely determined by them. These solitons are diverged with the velocity lni 2 '^) . 
They have varying amplitude and width in the general case, but we present also exam- 
ples where amplitude and width are constant. In these cases curved and weakly curved 
asymptotic solitons are both constructed. The lines of constant phase of the weakly curved 
solitons are deviated from the straight line just on a value proportional to lny 2 . 

Asymptotic solitons ( |1.5|) of the JE-I and the KP-I (j|], fll|| ) coincide taking into 
account transformations (|1.3|) and (FT 



2 Existence of a Johnson equation solution. 

After application of the change of variables 

to the scheme of the V.E.Zakharov and A.B.Shabat "dressing method" |M| of integration 
of the KP we obtain the following formulas for the JE. A JE solution has the form 

v(x,y,t)=2^-K(x,x,y,t), (2.1) 
where the function K(x, s, y, t) is a solution of the Marchenko integral equation 

poo 

K(x,z,y,t)+F(x,z,y,t)+ K(x,Z,y,t)F(£, z,y,t)<% = 0. (2.2) 

J X 

This equation is an equation with respect to z and x,y,t are parameters. The kernel 
F(x, z,y,t) of (|2.2j) satisfies the system of linear differential equations 



2 

h + 4g Q 2 ^ x + + ~^^ xx ~ ^ zz ^ ^ xxx ^ zzz = ® 



aFy + ^ {Fx + Fz) + I {Fxx ~ Fz ' 
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(a = i for the JE-I, and a = 1 for the JE-II). The correspondence of (|2.l| )- (|2.3D to the 



JE can be verified by a direct method, described in [32| for the KP. The scheme ( |2.1| )- 
( |2.3| ) don't allow us to solve the Cauchy problem for JE, but it is rather convenient for 
the construction of classes of solutions with various properties, in particular, of rapidly 
decaying rational solutions described in B and of non-decaying solutions. 
A wide class of solutions of (^|) as a = i can be found as follows 

F(x,z,y,t)= exp[ip(x- z)- q(x + z)+2qf{p,q,y)t]dfi{p,q), (2.4) 
J J n 

where Q, C C + (C + = {A | A = p + iq, q > 0} is the upper half-plane of the complex plane), 
f{Pi y) = Q 2 ~ 3p 2 + py/2 — y 2 /48, and dpi(p, q) is some measure on fl. 

To construct a JE-I solution by the scheme (|2.l| )-( |2~4|) we must define the set f2 in ( |2.4[) 
and the measure dfj,(p, q) on this set. For this goal we introduce two positive functions C(s) 
and g(s) which play an important role in the construction of the solution and investigation 
of its asymptotic behaviour, and we formulate the following conditions. 

Condition A. The function C(s) : R — >• M + is of class C 2 and such that 

C{s)>S>e 2 (<$,£ = const > 0), C"(s) > -1/24. (2.5) 

Condition B. The set £1 has the form 

n = {{p,q) GM 2 |-oo<p<oo, 0<£<g< h{p)} , (2.6) 

where q = h(p) is the envelope of the family of hyperbolas 

f(p,q,s) = C(s), (2.7) 

which touch it at the point 



(po(y),qo(y)) = ( C'(y) +y/12, y/c(y) + 3 (C'(y)Y 



Remark. The special structure of (|2.6| ) implies that 

C(s)= max f(p,q,8). (2.8) 

Condition C. The function g(s) : R — > M + , g(s) < A = const is of class C°° and 
such that the measure d\x of the form dfi(p,q) = g(p, q)dpdq with real positive g G C°° , 
g(po(Y), qo(Y)) = g(Y), satisfies the inequality 

Va = const > : // e a( - q+ \ pn dfi(p,q) < oo. (2.9) 
J J n 



Let us show that under Conditions A-C the scheme (|2.1|)-( |2.3| ) determines a smooth 
real solution of the JE-I vanishing as x — > +oo. For a function h(y) E L 2 [x,oo) define the 
operator F by 

POO 

[Fh](z)= F(s,z,y,t)h(s)ds (2.10) 



with the kernel F given by (2.4), where h(s) also depends on the parameters y,t. 
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Lemma 1. Assume that Conditions A-C are fulfilled. 

Then F is a self-adjoint, compact and positive operator in L 2 [x,oo). 



Proof. Its self-adjointness follows from the form ( |2.4| ) of F(x, z,y,t). Let us show that 
this operator is compact. We estimate the Hilbert- Schmidt norm of F: 



\F\ 



L 2 [x,oo) 



< 



< 



\F(x, z, y, t)\ 2 dsdz 
^2 J J dviP'i) J J ex P[ 4, 7(W + f(P,q,y))t]dfj,(p,q) 



4e 2 
1 

4^2 



d(J,(p, q) 



dfj,(p, q) < oo. 



Thus F is a Hilbert-Schmidt operator. Hence, it is a compact operator ( ||33|1 ). To prove 
the positivity of F, consider the scalar product 



(Fh, h) 



F(s, z, y, t)h(s)dsh(z)dz 



3 2<?/(p,(?,y)t 



,{ip-q)s 



h{s)ds 



dn(p, q) > 



as h(s) ^ 0. 

Under the conditions of Lemma 1 the following statement holds. 



Lemma 2. The scheme (2.1 )-(2.4) determines a smooth real solution of the JE- 1 vanishing 
as x — > oo and bounded for all fixed x, y,t (t > 0). 



Proof. Let us represent (2.3) in the operator form in L 2 [x,oo) 
ip + F<p = f, 



(2.11) 



where F has the form ( [2.10 ), and ip = K(x, z,y,t), f = —F(x, z, y, t). Due to the positivity 
of F the homogeneous equation cp + F(p = has only the trivial solution. Since F is a 
compact operator, then by the Fredholm theorem ([33]) inhomogeneous equation ( p. 11 ) 
has a unique solution given by 

K(x, z, y, t) = -(I + F)^ 1 F(x, z, y, t) 

with IK/ + F)- 1 ]] < 1. 

Due to Condition C and the fact that Q is inside the upper half-plane at positive 
distance from the g-axis (Condition B), F is an infinitely differentiable function with 



respect to all variables. Moreover, Df F — > (D c 



a 



0, 1, . . . ; i = 1, . . . , 4) as 



x + z — > oo, and DfF are bounded for all fixed x, z, y,t (t > 0). One can show (]|34|[) that 
the function K has the same properties. 

Let us prove that K(x, x,y,t) is a real function. After multiplication of ( 2.1 1) ) by Tp 
and integration with respect to z from x to +oo we obtain 



\\<p\\ 2 + (Fip,p) = (f,p). 



(2.12) 
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The self-adjointness of F implies that the imaginary part of the left-hand-side of ( [2.12 ) is 
equal to zero: 

Im / F(x,s)K(x,s)ds = 0. (2.13) 

J X 

Application of conjugation to ( |2.3|) for z = x gives 

POO 

K(x,x) + F(x,x) + / K(x,g)F(x,(;)dt = 0. 

J X 

It follows from ( p. 13 ) and the reality of F(x, x, y, t) that K(x, x, y, t) is real. ■ 



3 Theorem about long-time asymptotic behaviour of JE-I 
non-decaying solutions 

Our goal is to investigate the long-time asymptotic behaviour of the JE-I solution defined 
in the previous section. To define a domain in which we shall carry out the investigations, 
we introduce the following definition. 

Definition. Let M > 2 be an arbitrary number. The domain C K 2 given by 

G M (t) = l(x,y)eR 2 | \lng(y)\ < lnt, x > C(y)t - —L^ lnt M+1 

is called the neighbourhood of the solution front. 

The following theorem describes the asymptotic behaviour of the JE-I solution defined 
by Lemma 2 for large time. 

Theorem 1. Assume that Conditions A-C are fulfilled. 

Then the JE-I solution v(x,y,t) constructed by the scheme ( |2.lD ~(^7^) is represented in 
the domain Gm{P) as t ^ oo in the following way 

[M-l] 

v(x,y,t) = Vn(x,y,t) + 0/- T ^— j , (0<ei<l/2) (3.1) 

n=l ^ ' 

v n {x,y,t) - 



cosh 2 



q (Y) (x - C(Y)t + — ^ (lntn+l/2 - \ng(Y) - ln0 n (y) 



where q (Y) = y/c(y) + 3 (C'(y)) 2 , 

(C(y) + 48(C'(y)) 2 ) n - l (l + 24C" \y)) n - l l 2 Q^T^ 



2 (2n+5)/2(( n _ l)!)2(C(y) + 12(C / (y)) 2 )( 10 ™- 3 )/ 4 Q("- 1 )r(™- 1 ) ' 
and r( n ),Q( n ) > are the determinants of the n by n matrices with entries 

rfiUi = r ( 1± | ±i ) C 1 + = ^ + k + 1), 

i, k = 0, . . . , n — 1. 

Here the asymptotic representation Q3.1| ) is uniform with respect to x and y in Gm(4) for 
any fixed M > 2. 
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Let us mark the key points of the proof. The proof consists in three steps. On the 
first step we show that as t — > oo the kernel F(x, z,y,t) of integral equation ( |2.2| ) is 
represented as the sum of a degenerate kernel and a kernel with small operator norm in 
the space L 2 [x, oo). On the second step we prove that the degenerate kernel gives the main 
contribution in the asymptotic representation of the solution of equation ( |2.2| ). The third 
step consists in the analysis of representation ( fO| ) for the solution v(x,y,t) as t — > oo, 
where the function K(x, z, y, t) is a solution of the Marchenko integral equation with the 
degenerate kernel. 



First step. To investigate the Marchenko equation kernel F(x, z,y,t) (|2.4| ) as t — > oo, we 
set x = C{y)t + e, _z = C(y)t + C and define F(£, C, y, t) = F(£ + C(y)t, C + C(j/)t, y, i). 
Then the function F(£, £, y, t) is written as follows: 

Ht, C, y, t) = // ex P [^(e - c) - <z(£ + C) - 2<?(C(y) - f(p, q, y))tW(p, q). 

For sufficiently small e' > let us consider the curve 

2g(/(p,<z,y)-C(y))+ £ / = 0. (3.2) 

This curve separates the domain £1 into two subdomains O e > and Vt £ i so that f2 = O e /Uf2 e /. 
Here O e / lies between the curve q = h(p) and curve (|3.2|), moreover (po(y), Qo(y)) S 
O e '. The set ri e / is the complement of the set O e / in the domain O. According to this 



decomposition, kernel (^4) is the sum of two kernels which we denote F\(x,z,y,t) and 
F2(x, z, y, t) respectively. 

Let us make a change of variables, setting 

r = 2q(C(y)-f(p,q,y)) (3.3) 

in the kernel F\{x, z, y, t) which contains integration over the set O e >. Let u be the pro- 
jection of a radius vector directed from the point (po(y), <Zo(y)) to the point (p, q) G O e > 
on the tangent to the curve h(p,q) = at the point (po(y),qo(y)) or, that is the same, on 
the tangent to the curve f(p,q,y) = C(y) at the same point, i.e.: 

I2p (y) -y ( , y. 4^ (y) . , ^ 

u = — (g — Qo[y)) H — (p — Po(y))- 

Vi6^(y) + (i2 P0 (y)-y) 2W V^y) + (i2po(y) - y) 2 

(3.4) 

The system of equations ( |3.3| ), ( |3.4| ) has a unique solution with respect to p and g in 
O e as e' < 2 ^=. Therefore in the neighbourhood of the point (po>?o) the variables p and g 
can be expressed via the variables r and u: 



p(r, u) = po + fcxr + + k^ur + ^r 2 + k$u 2 . . . , 
g(r, u) = qo + Air + A2U + \?,ur + A/ir 2 + A5U 2 . . . , 
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where k n , X n are the coefficients of the corresponding Taylor series. The first have the 
form 



h(y) 

Ai(y) 
A 2 (y) 



i2po(y) - y 



( ?0 (y)(16g 2 (y) + (12p (y)-y)2)' 
4go(y) 

^16g 2 (y) + (12;po(y)-y) 2 ' 
4 

~16g2( y ) + (i2po(y)-y) 2 ' 
lggoCy) - y 

Vl6g 2 (y) + (12p (y)-y) 2 ' 



One can obtain the expansion coefficients k n and A n in an explicit form after n-times 
differentiation of ( |3.3j ) and fl3.4| ) with respect to r and u. In the neighbourhood of the 
point (po(y), Qo(y)) the equation q = h(p) can be written using variables r and u. It is 
easy to check that ^h(p)\ p=po ^ / since C" > —1/24 (Condition A). Therefore the 
curves f(p, q, y) = C(y) and q = h(p) have a contact of the first order, and q = h(p) takes 
the form u = u(r): 

u = ±a(y)y/r + b(y)r + . . . , 

where 



a{y) 



h pp (po 



16g 2 (y) + (12po(y) - y) s 



2?o(y) (48<? 2 (y) - (12po(y) - y) 2 - 16V(Po(y))^(y)) 

d 2 %) . 



dp 



2 \P=P0- 



Using the new variables u and r and the notation Eq(£, £,y) = e vpo ^^ Q 9o«/)(£+C) we 
write the function Fi(£, (, y, t) = F\{^ + C(y)t, ( + C(y)t, y, t) as follows: 

Fi(t,(,y,t) = E Q (£,(,y) 

x r dr r r+ "' dni(r,M,y)y(r,n,y)e i(p - po)( ^°- ( ' ? -' 2o)( ^ ) " ri (3.5) 

io J—a«/r+... 

where j(r,u,y) = j(p(r,u,y),q(r,u,y)) is the Jacobian corresponding to the change of 



variables (p, q) — ► (r, u). Let us expand integrand in (3.5) into a series with respect to the 



powers of r and u in the neighbourhood of the point (po(y), Qo(y)) (u = 0, r = 0): 
j( r >^y)5(^ w,y)exp[i(p-p (y))(£ - C) - 0? - ?o(z/))(£ + C)] 

oo n j n— j 

= EEE E C i r- J r i+m ^- / - m ^, m (y)(l + ^ n (r,n)), 

n=0 j=0 «=0 m=0 
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where 



<Pn,j,i, m {y) = 77-77 M — \ i / ■ n . {iki(v) + \i(y)) l (ik 2 (y) + Hy)) 1 " 1 

L\m\[n — J — m)\(j — 1)1 



Hy)=j(0,0,y) 



x(ih(y) - X 1 (y)r(\ 2 (y) - ik 2 (y)r~^ m g(y)j (y), 
1 



qo (y) ViGg'^y) + (i2 Po (y) - y) 2 ' 



g(y) = g(0,0,y), and the functions ifi n (r,u) satisfy \ip n (r,u)\ < An(r + \u\). 

After integration with respect to u and r, and performance of natural estimates, we 
obtain 

JV-lJV-n-1 , / \ 

F 1 (t,c,y,t) = E ti,cy)J2 E c n e j t( r+i+3)/ 2 o- + **(*)) + A ^ C, v, t), 

n=0 j=0 

where 

^ w = gMMgg^M r ( 2±i±l) (i + (_!)»«,, p.e) 



and |(5 n (t)| < 



2n\j\ 



'IKS 



\A N \ < A(N)g(y)J2 ^j^^KM , (3.7) 

3=0 

( e N+3 \ 1/2 
where A(N) < ( — - 3)jV+32jV+4 j . Estimate Q is valid as |£| < t 1/4 , |C| < * 1/4 - 

Let us estimate now the second kernel -^2(£>C>2/>*) containing the integration over the 
set £ /. Taking into account that C(y) — f(p,q,y) > §^ ((p,q) G £Vj <? > £ > 0) and 
Condition C, we can write 



*2(£,C,y,t)l< fL e-rt«W°M-f<r™W dii(p,q) 

J J il £ r 



- e ~^jJn e9m+Kl) d ^l) =0(e"^*) • (3-8) 

Thus assuming |£| < i 1 / 4 , |£| < i 1 / 4 we obtain the final asymptotic formula 

N-lN-n-l , v 

F(x,z, y ,t) = £;o(e,C,y)E E Jffffi/a (! + *»(*)) 

ra=0 j=0 

+ o(EP3^ M4+C, )+o( e -*') (3.9) 

with {; = x — C(y)t, C = x — C(y)t and |<5 n (i)| < i~ 1//2 . It is not difficult to see that this 
asymptotic expression can be differentiated with respect to x. 
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The following estimates hold: 

A 2 (AO 



f°° f°° A (N) 

J J lA^z^^tfdsdz^-fj^, (0< ei <l/2) (3.10) 



' X J X 

in the domain 



1 t M 2N + 5 

2go(y) y(y) 4 

and 

poo roo 

/ / \F 2 (s,z,y,t)\ 2 dsdz = 0(e- e,t ), (3.11) 

in the domain x > C(y)t — \ft. 

Taking into account (|3.9|) , (|3.1C| ) and ( 3.11| ) we can formulate the following lemma. 

Lemma 3. Inside the domain 

C>£>-^ylnt M+ \ \lng(y)\ < lnt, M > 1, (3.12) 

as t — > oo i/ie kernel F(£, £, y, £) is represented in the form 
F(£, C, y, t) = F N (£, C, y, t) + G(£, C, y, t), 

w/iere 

N-lN-n-l , , v 

^(e,c,»,*) = e*»tt-o-*tt + «x; E cTtS^ 



t (n+j+3)/2 ' 
n=0 j=0 

A/" = [(4M — 5)/2], and the functions ip n j{y) are bounded and defined in ( |3.6| ). The function 
G(s, z,y,t) admits the uniform estimate with respect to y in ( 3.12j ): 

/■OO /"OO / 1 \ 

J J \G( S ,z,y,t)\ 2 dsdz = Oi^ Y j^ ; j (0< ei <l/2). (3.13) 

Second step. We show that after replacing the kernel F(x, z, y, t) by the degenerate kernel 
Fn(x, z,y,t) one can obtain an asymptotic representation of the solution K(x,z,y,t) of 
the Marchenko integral equation ( fE2"| ) up to 0(t _1 / 2+£1 ) (0 < £\ < 1/2) as t — ► oo. Set 

up £ = x — C(y)t, C, = z — C(y)t and consider the domain z > x > C(y) -— lni M+1 , 

2<?o(y) 

| In g(y)| < lnt with an arbitrary number M > 1. Let us introduce the operators 

/*oo /'OO 

(F N f){z)= F N (s,z,y,t)f(s)ds, {G N f){z)= G N (s, z,y,t)f(s)ds, 

J X J X 

in L 2 [x,oo). Here F^{x, z,y,t) is the degenerate kernel 



N-lN-n-l , , \ 

^«,C,»,*) = ^° (€ - -* (€+0 E E ^ .(Ig^ , AT=[(4M-5)/2], 

n=0 j=0 



(3.14) 
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with £ = x — C(y)t, £ = z — C(y)t, and ip n j{y) defined by (|3.6f) . Gn(x, z,y,t) is the 
difference between -F(x, z, y, f) ( |2.4[ ) and F^(x, z, y, t): 

G N = F- F N . 

Now, ( |2.2| ) acquires the form 

(I + F N )f + G N f = h N + g N , (3.15) 

where / = K(x,z,y,t), h N = -F N (x, z,y,t), g N = -G N (x, z,y,t). By virtue of fl3.13|) 
one can easily obtain the following estimates for the norms of the operator Gm and the 

vector <7at in L 2 [x,oo) as z > x > C(y) t— rlni +1 , |lny(y)| < lni, i — > oo: 

2g (y) 

IIGjvH < ^it" 1/2+£1 , ||5jv|| < iir 1/2+£l (A 1 = const, < e 1 < 1/2). (3.16) 

The operator / + Fn is the direct sum of the two operators I\ + Fn and The first one 
acts in the subspace H\ of L 2 [x, oo), which is generated by the vectors 



e 



(ip -q )z ze (ipo-qo)z z N-l e {ip -q )z 



The second operator I 2 (I = I\@ I2) acts in the orthogonal complement H 2 = L 2 [x, 00) 
H\. Since we have F/v = F — Gn and the operator / + F is invertible, we deduce that the 
operator / + Fn is also invertible in L 2 [x, 00) and 

+ Fn)- 1 ]]^^ < A 2 . (3.17) 



We shall look for a solution of ( 3 . 1 5| ) of the form / = 4>n + tpN, where 4>n is the solution 



of the equation (/ + Fn)<Pn = hN. It implies that ipN satisfies 

(I + F N )ip N = g N - G N <fi N . 
According to (|3.17[) , we have 

I foMI <M\\9n\\ + ||Giv|||hMI). 

It follows from considerations presented below that 4>n is uniformly bounded with respect 
to (x, y) S Gjif(t) and t in the space L 2 [x, 00) n C[x, 00). This fact and ( 3.16| ) allow us to 
conclude that as t — > 00 

f(z) = ( f )N (z) + 0(t- 1 ^), 0< £l <l/2. (3.18) 
Third step. The replacement of the kernel F by the degenerate kernel F^ in equation ( |2.2[) 



and the implementation of the substitutions x = C(y)t + £, z = C(y)t + £ allow us to 
obtain the following integral equation for the function 1£jv(£j C> 2/> = + C(y)t,Q + 
C(y)t,y,t) (OO- 

i^v(£, a, y, i)Fjv(s, C, y, t)ds = 0, (3.19) 
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where Fn is given by ( 3.14| ). 

According to (2.1) and ( |3.18| ), the following representation of the function v(x,y,t) is 
valid in Gjvf(t): 



d 

v[x, y, t) = 2—K N {£, £, y, t) 



+ o(t- 1 / 2+£i ) 



(3.20) 



Z=x-C(y)t 

We shall look for a solution of equation (|3.19| ) of the form 
N-l 



-(*Po+<?o)C 



(3.21) 



n=0 



After substitution of flOID into (jjTIsD we obtain a system of algebraic equations for the 
function J n (^,y,t): 



N-l N-n-1 , , \ 

In f l"i ^ +(n+j+3)/2 

m=0 j=0 



s j+m e- 2goS ds 



N-n-1 

^ +(n+j+3)/2^ C 



3=0 



(3.22) 



n = 0, ...,JV-1. 

The solution of ( 3.22] ) has the form 

_ det[I + A(H,y,t)p 
11 det[I + A(S,y,t)] ' 

where I is the identity matrix, A(£,y,t) is the matrix with entries 

N-n-1 



[A] 



j=0 



nj 



(3.23) 



L 3+m 



s j+m e- 2qos ds, n,m = 0, . . . ,N — 1. 



The matrix [I + A(£,y,t)]( 1 ' is obtained via the substitution of the column of right-hand 
sides of the system ( |3.22| ) instead of Z-th column of the matrix [I+A(£, y, t)]. The functions 
ip n j are defined by (|3.6| ). 

The substitution of 7 n (£, y, t) into (|3.2l| ) gives us 



JV-l 



v tt /- *\ ST det(/ + A)W 

K N (U,y,t) = j: Ce 

n=0 



n„-(ipo+qoK 



Hence for the solution v(x,y,t) the following asymptotic representation holds true 

d 2 i 
v(x,y,t) =2-^\ndet[I + A(x-C(y)t,y,t)]+0(t-2 +£l ) (t^oo). (3.24) 
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To obtain asymptotics (|3.24|) we need to investigate the behaviour of A = A(£,y,t) = 
det[J + A(£,y,t)) in the domain as t — > oo. An analysis of the structure of the 

matrix / + A(^, y, t) shows that A(£, y, t) can be represented in the form 



N 



^y^ = ^E^Me- 2nq °^, (3.25) 



n=l 



£n(n+2)/2 



where P n (^,y,t) are polynomials with respect to £ of degree at most N n with coefficients 
bounded with respect to y and t. We have the following asymptotic relations as n < 
[(N + l)/2] 

^nl2)/2 e ~ 2nq ° m = det[AM(6y,t)](l + 0(r^)), (3.26) 

where A^ n \^,y,t) are the n by n matrices with entries 
n— 1 , j 

(n) _ \ -» Yij-Lj+k 



A {n > ~ V 



^(i+J+3)/2 ■ 
j=0 

The matrix A' n '(^, y, f) is obviously written as the product of two matrices. Therefore we 
have from (|3.6|) and ( |3.23 ): 

**[A«>&y,t)] = £m®!!?to ^ n)iH + fy\ (3.27) 
where ftW and / (n) (£,2/) 

are the determinants of the n by n matrices with entries 

r oo (3.28) 
J (n) - 1-^ - / s i+fc e _2s ds 
J io{y)£, 

i, k = 0, 1, . . . , n — 1. 
The elements of the determinant /( n ) satisfy the relations 

^h±^- = -2q (y)I i+k + + k)q Q (y)Ii+k-i, 
whence we have the equality 



n— 1 



+ 2ng (y)/ {n) =go(y)E i? fc ) i' 



fc=0 



where -B^."\ is the determinant obtained from 1^ by replacing the (k + l)-th column by 
the column with entries ili+t-i (i = 0, 1, . . . , n — 1). We can show that X^fc=o ^k+i = 
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hence I {n) (£,y) = l^{0,y)e- 2nqo ^. It follows from (|^28|) that l {n) (0,y) do not depend 
on y, therefore we write finally 

where Iq are positive numbers depending on n. 

The determinant Tq 7 ^ is, up to a factor 2~" 2 , the Gram determinant of the system of 
functions x k e~ x / 2 (k = 0, 1, ... ,n — 1) on the semi-axis [0,oo), therefore / 0. The 
determinants f^ n ) are also the Gram determinants of the system of functions 



' (_ 1 )(fc+l)/4 a .(2fc-l)/4 ea! /2 j _oo < x < 
a .(2fc-l)/4 e -x/2 < X < OO 



on the axis (— oo, oo), and fiH / 0. Since in the region ( |3.12|) the determinant A does not 
vanish (the operator (J + F)' 1 is bounded), it is not difficult to see, using ( 3.25| )-( |3T27| ) , 
that fl^I^ > 0. 

Thus, for n < [(N + l)/2], P n (^,y,t) = P n (y,t) does not depend on £, where 



P 



{y, ) ~ 2^ qo{ y)r^ k -j o{k ir {1+G{t )} 



> 0. 



(3.29) 



Let us turn back to the representation (|3.24| ) of the JE-I solution. We can write 

d 2 A" A- (A') 2 
v(x,y,t) ~ 2^ In A = 2 ^J-. 



(3.30) 



Let us cover the domain Gjv/(t) (M > 2) by the subdomains 



a i 



1 t 2+e 
-- — In— r^- < £ < oo \ , 



2go ff(y) 



l , t n+1+e l , t n - £ 

■ In — < £ < - - — In 



2<?o g(y) 



2<?o g(y) 



n = 2, 3, . . . , m — 1, 



In —r^ < (, < — - — In 



where m 



2<?o 5(2/) 



2<?o 5(2/) J ' 



[M — 1], £ = x — C{y)t. With £ being inside any specific a n the 



numerator and denominator in ( 3.30| ) have the asymptotic representations 

^ J f (n-l)(n+l)/2 i(2n+l)/2 ^ ^ ^ 

2 



A 



; ff 2(n-l) e -4(n-l)g (?/)£ 
f(n-l)(n+l) 



-29o(t0£ 



t(2n+l)/2 



;i + o(t- 1 /2)) 



(3.31) 
(3.32) 



with 



^n(y) 



jS(y)a n2 (v)rt n) lo n) 

2"(g (y))" 2 re=o(^o : 
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Using (jOH) and (|3T32b , we conclude that the solution of the JE-I possesses the following 
asymptotic representation uniformly with respect to £: 



v(x,y,t) 



8q (y)gP n -l-?n t -(2n+l)/2 



P _l_ r,P e" 2 90(»)g 

* n— 1 ~r i/^n £-(2n+i)/2 



> J— t (2n+l)/2 



1 + 



P„_l t(2" + l)/2 



£=*-C(y)t 



2% 2 (y) 



cosh 



(y) (x - C(y)t + (ln*n+V2 - ln 5 (y) - In 



*~»v»/ *~p B _i(y), 
Thus, the JE-I solution splits in Gj^it) into [M — 1] solitons of the form 

2g 2 (y) 

v n (x,y,t) = — 

cosh [<7o(y)<Mz,y,t)J 

with curved lines of constant phase 

Mx,V,t) = x - C(y)t + — (mt n+1 / 2 - In g(y) - ln0 n (y)) 

where 

jo(y)« 2n - 1 (y) 



0n(l/) 



2(go(y)) 2n - 1 [("-i)!] 2 ^-i)4 n - 1) 

(C(y) + 48(C / (y)) 2 )"~ 1 (l + 24g // (j/))"- 1 / 2 QWr(") 
2 (2n+5)/2(( n _ i)!)2(c(y) + 12(C / (y)) 2 )( 10 ™- 3 )/ 4 Q( n - 1 )r(™- 1 ) 



and Q( n \ Jq are the determinants of the n by n matrices with entries 



[n (n Ww = r f ) ( i + (-1^ 



[4 n) ]i+i,fc+i = T(i + k + 1), i, fc = 0, 1, . . . , n - 1. 
The theorem is proved. 



(3.33) 



(3.34) 



4 Examples 

Example 1. Assume that 



Then 



Po(2/) = g, 90(1/) = — ^ — • 
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and the set f2 has the form 



9. 



f / \ a/36p 2 + 1 
= |(Pj <Z) — oo<p<oo, 0<e<g< 

Determine y(y) and d/x as follows 

g(y) = e~y\ dfi = e^ 2 ^- 1 ' 2 ^ dpdq. 

Then Conditions A-C are fulfilled, and according to Theorem 1 there exists a JE-I solution 
which splits as t — ► oo in the domain (M > 2) 



G M {t) = {(x,y) 



|y| < vlni, x > — + — 



— lr 1 ^+ 1 | 

24 i6 y^TT J 



into [M — 1] curved solitons of the form 

y 2 + i 



v n (x,y,t) 



8 cosh 2 



^±^„(*,y,i) 



ipn(x,y,t) = x - 



y 2 t 



+ 



24 16 V V+I 



lnt "+l/2 +y 2 _ ln 



3 n - 1 / 2 (6y 2 + l)«-iQ(")r( n ) 



2 5 ™-3/ 2 (2y 2 + l)( 10 ™- 3 )/ 4 [(n - l)!] 2 Q( n - 1 )r( n - 1 ) 



Example 2. Suppose now that 

C(y) = b 2 (b = const > e > 0). 

Then 

Po(y) = ?o(y) = 

and the set £1 has the form 

^ = {(P; q) I — 00 < P < °°, < e < g < 6} . 
Determine y(y) and d/^ as follows: 



g{y) 



d\i = e 



dpdq. 



Then Conditions A-C are fulfilled, and according to Theorem 1 there exists a JE-I solution 
which splits as t — ► oo in the domain (M > 2) 



G M (t) = {(x,y)eR 

into [M — 1] curved solitons of the form 

2b 2 

v n (x,y,t) = — , 

cosh [bip n {x,y,t)\ 

ip n (x,y,t) = x-b 2 t 

+ ±(lnt n+1 / 2 +y 2 -ln 



1 



< v'ln/. , >6 2 t-^lnt M + 1 } 



Q(n) r (n) 



j3(n-l/2)[( n _ 1 )!]2Q(n-l)r(n-l) I ' 
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Example 3. We consider the function C(y) (and also the functions po(y), qo(y) and the 
set fi) introduced in the previous example. Such structure of Q allows us to formulate a 
weaker condition on the function g(y). We suppose that g(y) = (1 + y 2 ")^ 1 with integer 
a > 4 and the measure dfj, instead of ( |2,9| ) satisfies the following inequality: 



dn(p, q) ^ 

< oo. 



1 + (I2p) 2a 



Then the function v(x,y,t) constructed by the scheme (2.1)-(2.3) satisfies the JE-I, but it 
is not infinitely differentiable. We can prove that such a JE-I solution splits as t — > oo in 
the domain (M > 2) 

G M (t) = {(x,y)eM. 2 | y 2a <t, x > b 2 t - ^lnt M+1 ] 

into [M — 1] curved solitons of the form 

2b 2 

V n (X,y,t) 




cosh 2 [bip n (x,y,t)] ' 

ij) n (x,y,t) = 

In t n +V 2 + lnf 1 + v 2a ) - In " 

+ m { i+y ) ln fe3(n _ 1/2)[(n _ 1)!]2Q(n _ 1)r(ri _ 1) j • 

Their lines of constant phase are deviated from the straight line just on the value 2 hi(l + 
y 2a ). Therefore we call them weakly curved solitons. 

The approach developed in the present paper can be applied to solve inverse problems 
for other dispersion models, like the first-order Debye model, or more generally, the multi- 
resonance Lorentz and A^th-order Debye models. The generalization is straightforward 
and is based on the construction of proper piecewise holomorphic functions near each pole 
of the model. 
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